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ABSTRACT 

The  propagation  of  hydromagnetic  waves  is  examined 
In  a  perfectly  conducting  compressible  fluid.   The  medlvun 
is  highly  anisotropic  and  exhibits  strange  properties,  many 
of  which  are  hidden  by  Fourier  analysis  or  by  examination 
of  plane  waves.   In  a  certain  sense,  it  is  possible  to 
separate  an  arbitrary  wave  motion  into  three  modes  (part 
of  the  separation  is  always  possible,  the  remainder  only 
In  certain  limiting  cases).   However,  in  separating,  one  can 
lose  the  domain  of  dependence  properties  of  the  original 
hyperbolic  system.   One  of  the  modes  (in  certain  limiting 
cases,  two  modes)  represents  one-dimensional  wave  propaga- 
tion along  magnetic  field  lines  without  three-dimensional 
radial  attenuation.   Depending  on  the  amount  of  dissipation 
present,  such  signals  could  propagate  large  distances. 
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PROPAGATION  OF  MAGxNETOHYDRODYNAMIC  V/AVES 
V/ITHOUT  RADIAL  AT'.  EKUATI  Owi/ 
1 ,      Introduction 

We  shall  be  concerned  with  the  I'ollowlnp  magneto- 
hydrodynamlc  formulation: 


i£  +  div  (pu)  =  0 


(1)     ^   p  (i^  +  u-  \7  uj  +  a^  7p  =  i  curl  B  x  B 


curl(u  X  B),   div  B  =  0  , 

ist^/ 


This  system  is  essentially  that  of  Lundquist— '   and  is  the 
simplest  formulation  of  compressible  magnetohydrodynamics . 
It  follows  from  the  removal  of  all  dissipative  mechanisms; 
more  precisely,  the  fluid  is  isentropic  ivith  a  scalar  pressure, 
and  Ohm's  law  is  taken  in  the  simple  form   E  +  uxB   =0  , 
Also,  displacement  current  and  electrostatic  forces  have  been 

dropped;  thus  the  system  (l)  is  Galilean  invariant  and  of 

3/ 
ccnservatirn  I'orm„— '   V/e  are  using  rationalized  MKS  units 

This  report  contains  the  text  of  an  address  presented  at 
the  Third  Lockheed  Symposium  on  Magnetohydrodynamics  on 
November  21,  1958. 

The  author  wishes  to  express  his  appreciation  to  Jack  Bazer, 
K.  0.  Frledrichs,  and  Herbert  Kranzer  I'or  valuable  dis- 
cussions on  this  subject. 

2  _, 

3.  Lundquist,  Studies  in  magneto -hydradynamic s ,  Arklv  fur 

Fyslk,  Band  5 ,  ( 1952  j . 
3 

For  a  discussion  of  these  points  see  A.  A.  Blank  and 

Harold  Grad,  Fluid  Magnetic  Equations  -  General  Properties, 

NYO-6J4.86,  Notes  on  Magneto -Ply  drodynam.lcs  VI,  Inst,  of 

Math.  Sci.,  N.Y.U.,  (July  1,  1958). 
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for   B;   p   and   u   are  the  gas  density  and  flow  velocity 
vector,  and   a   =  ^p/{)p   is  the  ordinary  gas  sound  speed. 

The  theory  of  this  nonlinear  system  has  been  shown  to 
be  basically  similar  to  that  of  ordinary  compressible  gas 
dynamics.   In  particular,  one-dimensional  simple  wave  and 
shock  wave  solutions  have  been  iiivestigatedo—   Neverthe- 
less, in  detail,  the  ma gnetohydro dynamic  phenomena  are 
decidedly  more  complex,.   One  distinction  is  that  there  are 
three  characteristic  (or  propagation)  speeds  for  the  magneto- 
hydrodynamic  system  in  contrast  with  the  single  sound  speed, 
a,   of  gas  dynamics.   What  is  probably  more  significant, 
however,  is  the  strongly  anisotropic  natiore  of  the  electri- 
cally conducting  gas.   In  fact,  two  of  the  tliree  propagation 
speeds  reduce  to  zero  in  a  direction  perpendicular  to  the 
magnetic  fields   This  is  sufficiently  degenerate  to  render 
inapplicable  most  of  the  conventional  theory  of  wave  propa- 
gation in  anisotropic  media. 

Most  of  the  analysis  will  refer  to  the  linearization  of 
(1): 

K.  0.  Priedrichs  and  Ho  Kranzer,  Nonlinear  Wave  Motion, 
NYO-6i4.86,  Notes  on  Magneto -Hydrodynamics  VIII,  Inst,  of 
Math.  Sci.,  N.Y.Uc,  (July  31,  19i?8), 

Jack  Bazer,  Resolution  of  an  Initial  Shear  Flow  Disconti- 
nuity in  One-Dimensional-Hydromagne tic  Flow,  Research 
Report  No.  "HH-5,  Dlv,  of  Electromagnetic  Research,  Inst» 
of  Math.  Scl.,  N»YcU.  (June,  1956). 

W.  B.  Ericson  and  Jack  Bazer,  Hyaroraagne tic  Shocks, 
Research  Report  No.  MH-S,  Div.  of  Electromagnetic 
Research,  Inst,  of  Math.  Sci,,  N.Y.U,  (January,  1953), 
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4f  +  p„  dlv  u  =  0 


'2'    {   Po^*%  VP=i 


curl  B  X  B 

[1  o 

i|  =  curl  (u  X  B^),   div  B  =  0  . 

The  unperturbed  quantities  are  denoted  by   p  ,  B  ,  a  ,  u  =0, 
while   p,  u,  B  here  denote  the  perturbations^   B   repre- 
sents a  uniform  magnetic  fie  id. 

Since  this  system  of  equations  is  hyperbolic,  its  solu- 
tions can  be  described  in  terms  of  the  propagation  of  waveso 
As  in  the  case  of  the  ordinary  wave  equation,  one  can  investi- 
gate a  simplified  problem,  viz.  the  propagation  of  wave  fronts 
and,  associated  with  this,  one  can  construct  a  theory  of  ray 
optics,—'   An  important  point  to  realize  is  that  the  first 
order  ( Hamilton- Jacobi )  partial  differential  equation  which 
governs  the  behaviour  of  wave  fronts  does  not  describe  the 
entire  content  of  the  second  order  wave  equation^   In 
particular,  there  may  be  a  complicated  non-steady 
flow  left  behind  a  wave  front  propagating  into  still  air, 
and  this  residue  can  only  be  analyzed  by  the  use  of  the 
full  wave  equation.   The  situation  is  more  complicated  in 
the  magnetohydro dynamic  case  where  there  are,  in  general, 
three  distinct  wave  fronts  which  propagate  at  different  speeds, 

Jack  Bazer  and  0.  Fleischman,  Geometric  Hydromagnetics , 
Research  Report  Nc  MH-12,  Div,  of  Electromagnetic 
Research,  Inst,  Math,  Sci » ,  N.YoU.  (February  1959 )• 
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For  example,  a  wave  front  which  moves  at  one  of  the  two 
slower  speeds  can  leave  a  "residue"  which  propagates  ahead 
of  itself  as  well  as  behind. 

This  remark  introduces  the  major  points  of  investiga- 
tion of  this  paper:   first,  the  discovery  of  special 
solutions  which  propagate  at  a  single  one  of  the  character- 
istic speeds,  and  second,  the  question  of  the  resolution  of 
an  arbitrary  solution  of  the  linear  system  (2)  into  three 
"modes"  which  propagate  at  distinct  speeds.   In  one 
dimension,  a  complete  separation  can  be  obtained  using  the 
Riemann  invariants  of  the  linear  system  (2).   For  the  non- 
linear system,  in  one  dimension,  special  simple  wave 
solutions  which  propagate  at  a  single  speed  can  be  found. 
Separation  of  the  wave  solutions  of  the  linear  system  (2) 
can  be  only  partially  accomplished  in  more  than  one  dimen- 
sion, and  the  natvire  of  this  separation  is  quite  subtle 
(section  7)»   However,  it  is  easy  to  exhibit  a  large  class 
of  special  solutions  which  propagate  at  a  single  one  of  the 
characteristic  speeds  (sections  5  and  6). 

In  contrast  to  the  problem  of  the  separation  of  general 
solutions,  wave  fronts  can  always  be  separated,  and  even  in 
the  case  when  the  unperturbed  magnetic  field,   B  ,   is  not 


uniform. 


1/ 


See  footnote  on  page  6. 
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It  should  be  remarked  that  a  Fourier  expansion  of  the 
linear  system  (2)  immediately  yields  separable  modes „   This 
resolution  is  artificial,  howeverp  and  hides  some  of  the  most 
important  properties  of  the  solutions.   For  example,  a  Fourier 
resolution  does  not  easily  answer  the  question  whether  one  can 
create  an  initial  distrubance  which  is  confined  to  a  finite 
region  and  which  propagates  as  a  single  one  of  the  three  modes  „■=/ 
With  modes  defined  in  a  more  physically  appealing   sense  than 
as  Fourier  components,  the  modes  cannot  be  entirely  decoupled 
even  for  the  linearized  systemo 

The  fact  that  two  of  the  propagation  speeds  are  in  some 
sense  zero  in  a  direction  perpendicular  to  B    suggests  the 
possibility  of  waves  which  propagate  one-dimens ionally  with- 
out radial  spreading »   However,  the  situation  is  complicated 
by  the  interaction  between  different  modes  (eg»  a  v  i  ..  front 
of  one  type  can  leave  a  residue  involving  the  other  two  modes) 
and  by  some  ambiguity  in  the  meaning  of  "zero  propagation 
speed"  perpendicular  to   E    {associated  with  non-convexity 
of  one  of  the  characteristic  cones )o   The  analysis  made  here 
is  a  step  in  the  direction  of  mathematical  verification 
of  certain  qualitative  explanations  that  have  been  offered 
of  the  phenomenon  of  "whistlers",  but  there  is  considerably 
more  to  be  done  in  this  direction  (section  9)- 

A  Fourier  analysis  of  the  radiation  fromi  a  point  source  has 
been  made  by  Leo  C«  Levitt  in  his  thesis,  California  Insti- 
tute of  Technolog7/,  1957 »   However,  some  of  his  results  are 
difficult  to  interpret  because  of  domain  of  dependence 
ambiguities  o 
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2  o   Characteristics 

A  characteristic  surface   4)(x,y,z,t)  =constant,  (or  wave 
front)  of  the  system  (1)  can  be  defined  as  a  surface  in  space- 
time  across  which  there  can  exist  discontinuities  of  some  of 
the  first  derivatives  of   p,u,  or  B.   The  system  of  equations 
(1)  imposes  certain  relations  among  the  relative  magnitudes  of 
the  jumps  in  the  various  derivatives;  specifically,  it  reduces 
to  a  system  of  seven  (scalar)  homogeneous  linear  equations  in 
the  jumps  of  the  seven  derivatives  of  p,  u, B  normal  to  the 
surface^   Since  these  equations  are  homogeneous,  solutions 
exist  only  when  the  determinant  vanishes.   This  condition  re- 
stricts the  possible  orientation  which  can  be  taken  by  the 
surface  element  ^=   constant  at  a  given  point  (x,y,Zj,t)  and  is 
expressed  as  a  homogeneous  polynomial  in  the  derivatives 
(^  .'I'  ,^    J  4>4- )   which  define  the  normal  to  a  characteristic 
surface  element.   It  is  convenient  to  represent  this  normal  as 
a  vector  (n,  jn^  ,n-. , -c)   with  unit  projection  on  physical  space, 

pop  ~\/ 

nT'+nn+nT  =  lo      The    characteristic    equations    are    then—'' 
1       d.        J> 

(u    =-  c)6p    +   p5u„   =      0 
^   n  "^      n 

P  1 

(■:>)  I     p(u-c)5u+an5p  =  -(nx5B)xB 

1,     (u   -   c)5B  -u6B    -  B   5u  +B5u    =     0.      5B^   =   0      . 
\    ^   n  n        n  n  -  n 

We  have  written   5p   for  the  jump  in  the  normal  derivative  of 


^j_ 

Friedrichs  and  Kranzer,  footnote  on  page  5> 


9  - 


p  ,  etc.  since  one  can  show  that  the  same  characteristic  rela- 
tions hold  for  small  jumps  in  p  itself  as  for  finite  jumps  in 
^p/)n. 

It  should  be  remarked  that  the  system  (l)  is  mathematical- 
ly sound  even  when  the  restriction  divB=  0  is  not  imposed. 
A  property  of  this  extended  system  is  that  -vr(div  B)  =  0  for 
all  time.   This  reduces  the  system  (1)  to  sixth  order,  and  the 
characteristic  equations  (3)  to  a  set  of  six  independent 
equations . 

The  determinant  condition  is  a  homogeneous  polynomial  of 
the  sixth  degree  in   (n^ ,  n^,  n^,  -  c)  J 

ik)         [c^-  B^V^p][c^-(a2+  bVm.p)c^+  a^B^^/^p]  =  0,  B^=B.n 

The  linear  system  (2)  which  is  homogeneous  in  first  derivatives 
represents  a  non-dispersive  medium.   It  is  easy  to  verify  that 
this  characteristic  determinant  gives,  in  the  linear  case,  the 
speed  of  propagation,   c  ,  of  a  plane  wave  in  the  direction  n. 

It  is  convenient  to  represent  these  characteristic  speeds 
as  a  polar  plot  J  on  a  ray  with  direction  n.  ,  we  lay  off  radii 
equal  to  the  corresponding  values  of  c   ,  We  call  this  the 
normal  speed  locus.   For  the  ordinary  wave  equation,  it  is  a 
sphere  centered  at  the  origin.   For  the  present  equations,  the 
result  is  as  shown  in  Figure  1.   There  are  two  possibilities 
depending  on  whether  the  Alfven  speed, 

(5)  A   =   (B^  /np)^/2 
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A  <  a  A  >  a 

Figure   1 


A  =  a 

(b) 
Figure   2 


A  >  a 


(cj 
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is  larger  or  smaller  than  the  gas  sound  speed,  a  ,   The  trans- 
verse locus,  which  is  given  by  the  first  factor  of  (ij. )  ,  consists 
of  two  spheres  through  the  origin.   The  transverse  locus  sep- 
arates the  slow  locus  (which  also  has  two  lobes  5  meeting  at  the 
origin)  from  the  fast  locus  which  is  a  simple  convex  surface 
centered  about  the  origin.   In  any  given  direction,  n  ,  a  slow 
plane  wave  will  propagate  slower  than  the  transverse  wave  in 
the  same  direction,  and  this,  in  turn,  will  be  slower  than  the 
corresponding  fast  plane  wave.   For  A>a,  the  transverse  locus 
touches  the  fast  locus  when   n   is  parallel  to   B,  and  for  A<a, 
it  touches  the  slow  locus.   It  is  evident  that  plane  waves  of 
either  slow  or  transverse  type  propagate  with  zero  speed  when 
aimed  perpendicular  to   B„   Using  this  normal  speed  locus,  it 
is  possible  to  describe  a  simple  construction  (finite  differ- 
ence scheme)  to  obtain  a  characteristic  surface,   4i(x,y,  z,t) 
=  constant,  given  its  intersection  with  the  initial  manifold, 
t=0.   At  each  point  of  the  wave  front,  we  construct  the  normal 
and  lay  off  the  length  c  dt   where   c   corresponds  (on  the 
normal  speed  locus.  Pig,  1)  to  the  given  direction  of  the 
normal o   This  construction  merely  expresses  the  definition  of 
(n,  ,np,n,,-  c)   as  the  four-dimensional  normal  to  the  charac  = 
teristic  surfaces   The  construction  is  essentially  correct, 
but  as  we  shall  soon  see,  it  can  break  down  after  some  time. 
It  follows  immediately  from  this  construction  that,  if  an 
initial  wave  front  splits  into  slow,  transverse,  and  fast 
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components,  the  slow  front,  at  any  later  time,  will  not  have 
progressed  as  far  as  the  transverse  wave.   For,  at  the  initial 
instant  it  is  slower  than  the  transverse  wave  in  every  direc- 
tion, and  it  can  only  pass  the  transverse  wave  at  some  later 
time  by  touching  it  at  some  point,  and  at  this  presumed  point 
of  tangency  its  speed  would  again  be  less. 

For  many  purposes,  a  more  important  locus  is  the  one 
dual  to  the  normal  speed  locus.   Through  each  point  of  the 
normal  speed  locus  we  pass  a  plane  normal  to  nj  this  rep- 
resents a  possible  element  of  a  characteristic  surface 
4  =   const.   The  characteristic  locus  (Pig,  2)  is  the 
envelope  of  these  planes  5  ordinarily  one  speaks  of  the 
characteristic  cone  which  is  the  "surface"  in   (x,  y,  z^t) 
obtained  by  passing  a  ray  from  the  origin  to  each  point  of 
the  characteristic  locus  considered  to  be  drawn  in  the  three- 
space   t  =  l|  (the  term  "cone"  should  not  be  taken  too  liter- 
ally, since  the  envelope  can  be  quite  degenerate) » 

More  precisely,  one  considers  an  infinitesimal  charac- 
teristic cone  attached  to  each  point  (x,y,z,t)  ,  since  the 
characteristic  speeds  are  variable  in  general,   A  character- 
istic surface  has  the  property  that  it  is  everywhere  tangent 
to  the  local  characteristic  coneo   The  finite  characteristic 
surface,   (j)  =  cons tant ,  is  obtained  by  piecing  together  char- 
acteristic surface  elements  in  two  stages,  first  taking  an 
envelope  at  a  given  point  (x,y,  z,t),  then  taking  an  envelope 
of  these  cones, 
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For  the  special  case  of  linear  equations  with  constant 
co-efficients  (eg,  the  system  (2))^  the  characteristic  cone 
Is  Independent  of   x,y,z,   and   t.   Consequently,  the  evo- 
lution of  a  wave  front  from  time   t    to  time   t-,   can  be 
obtained  in  one  step  by  attaching  to  each  point  of  the 
initial  wave  front  a  characteristic  locus  whose  size  is 
given  by  the  scale  factor   t-,  -  t    and  then  taking  an 
envelope.   In  particular,  the  characteristic  locus  itself 
can  be  thought  of  as  the  wave  front  which  emerges  from  a 
point  disturbance  after  a  finite  time. 
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3o   Propagation  of  Wave  Fronts 

At  a  given  point   (x,y,z,t),   the  characteristic  system  (3) 
can  be  considered  to  be  a  set  of  linear  equations  with  c  as  an 
eigenvalue  parameter.   To  each  eigenvalue   c   (cf.  equation  {I4.)  ), 
there  exists  a  corresponding  eigenvector   (5p,6u,6B).   Expan- 
sion of  an  arbitrarily  given  set  of  jumps   (6p,5u,5B)   into 
eigenvectors  amounts  to  a  resolution  of  the  discontinuity  into 
six  modes  which  propagate  at  the  various  characteristic  speeds. 

Description  of  the  propagation  of  the  various  wave  fronts 
is  facilitated  by  introducing  the  rays  or  bi character! sties  of 
the  original  system.—   The  equation  for  a  characteristic 
surface  or  wave  front  is  (cf.  equation  (1^)): 

We  have  introduced  the  vector  Alfven  speed,   A  =  B/J\ip ,      for 
notational  simplicity.   This  is  a  first  order  partial  differ- 
ential equation  and  itself  has  characteristics  which,  because 
of  the  special  form  of  the  equation,  can  be  described  as  curves 
rather  than  as  cones.   For  simplicity  we  consider  (6)  to  be 
solved  for   c|).  =  F(7(|),a,A)  (six  roots).  The  characteristics  of  (6) 


^  =   <^F  ^   .^F  ,^a  ^  ^F  ^  U^ 
dt     ~  Jt        "^Ja  Jt    7a   Jt^ 


A  detailed  analysis  of  this  theory  is  given  by  Bazer  and 


Fleischman,  see  footnote  on  page  6. 
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These  equations  yield  not  only  the  rays   x(t)   but  also  the 
propagation  along  the  rays  of  the  characteristic  surface 
element   (}),  (t),  '74*  ^ ''' ^ '  ^^   the  special  linear  case  with 
constant  coefficients  ( i „ e c  (2)),  we  see  that  the  rays  are 
straight  lines  and  the  surface  element  is  carried  parallel 
along  a  rayo   In  the  characteristic  locus,  Fig.  2,  the  ray 
corresponding  to  a  given  surl'ace  element  is  the  line  from 
the  origin  to  the  point  on  the  locus  which  has  the  given 
tangent  direction.   The  image  of  a  wave  front  at  a  finite 
later  time  can  be  constructed  without  taking  envelopes 
merely  b"'  laying  ofl'  at  each  point  the  ray  of  finite  length 
(taken  from  Fig.  2)  which  corresponds  to  the  direction  of 
the  wave  front  at  the  given  pointj  in  addition,  we  know  that 
the  tangent  plane  will  be  the  same  at  the  end  of  the  ray. 
This  construction  should  be  compared  with  the  previous  one 
using  normals  to  the  front  and  the  normal  speed  locus ^ 

The  most  significant  portion  of  the  characteristic 
locus  (Fig.  2)  from  the  point  of  view  of  general  theory  is  the 
fast  locus  since  this  determines  the  domain  ol'  dependence.   This 
locus  is  a  simple  convex  surface  except  in  the  special  case 
when  the  Alfven  and  gas  speeds  are  equal.   In  this  case,  the 
envelope  proper  is  the  "outer"  part  of  the  apple-shaped  locus 
shown  in  Fig.  2b  (the  fast  and  slow  loci  merge  analytically 
in  this  special  case,  producing  the  entire  apple).   The 
vertical  dashed  line  (representing  a  disc)  is  a  singular  part 
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of  the  envelope.   It  is  required  by  the  general  theory  to  make 
the  locus  convex,  and  it  is  also  the  correct  limiting  locus  as 
A  — — >  a.   However,  this  part  of  the  characteristic  locus  is 
characterized  by  a  two  parameter  conical  family  of  rays j  there 
is  no  unique  ray  which  belongs  to  this  surface  element »   Although 
this  part  of  the  locus  is  a  carrier  of  discontinuities,  it  is 
not  a  carrier  of  the  dominant  discontinuity.   For  example,  the 
image  of  an  initial  wave  front  across  which  first  derivatives 
of   (p,u,B)   jump  will  have  continuous  first  derivatives  and 
will  have  a  higher  order  discontinuity  which  propagates  along 
these  rays  (Fig.  3)« 

The  transverse  locus  is  entirely  singular.   All  plane 
elements  pass  through  the  two  points   x=+  A   (the  B-direc- 
tion  is  taken  as  the   x-axls)<,   The  "dominant"  part  of  the 
locus  consists  of  these  two  points,  but  the  line  joining 
them  is  also  a  singular  envelope.   This  part  of  the  locus  is 
characterized  by  non-unique  rays  (  a  one  parameter  family). 
Using  the  transverse  characteristic  locus  to  obtain  the 
motion  of  a  wave  front,  one  obtains  the  apparently  one-dimen- 
sional undistorted  propagation  shown  in  Fig.  l^..      The  wave 
front  is  forever  contained  within  the  flux  tube  in  which  it 
originated o   The  dotted  part  of  the  construction  results 
fron  the  line  joining  x=+  A   in  the  characteristic  locus. 
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Figure   l\. 
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The  regular  part  of  the  slow  locus  is  not  convex 

and  has  the  two  cusped  lobes  shown  In  Fig.  2,   Again, 
there  is  a  singular  part  of  the  envelope  joining  the  two 
inner  cusps.   The  distance  of  one  of  the  outer  cusps  from 
the  origin  is  greater  than  the  maximum  distance  from  the 
origin  of  the  entire  slow  normal  speed  locus.   In  other 
words,  a  slow  disturbance  can  propagate  faster  than  the 
speed  of  any  slow  plane  wave,  and  it  can  also  propagate 
faster  in  some  directions  than  a  transverse  distxarbance. 


Figure  5 
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In  Fig»  5  we  sn.ow  the  evolution  of  an  orlf-inaliy 
spherical  slow  wave  front.   For  a  certain  time,  the  wave 
front  is  distorted  but  is  not  singular.   During  this  time, 
the  slow  wave  lags  behind  the  transverse  wave  and  remains 
within  the  original  flux  tube  (consistent  with  the  earlier 
construction  based  on  the  normal  speed  locus ) a   Howeverj 
cusps  develop  after  a  finite  time,  and  the  slew  front  can 
then  overtake  the  transverse  front o   In  particular,  the 
front  propagates  three-dimens ionally ,  penetrating  outside 
the  original  flux  tube. 

The  various  portions  of  the  characteristic  locus 
represent  poss  ible  carriers  of  disc ontinuities «   With 
special  initial  data,  part  of  the  locus  is  dropped,.   We 
shall  see  later  that  the  line  joining  the  transverse 
points   X  - +_   A   is  never  presento   It  has  not  been  veri- 
fied whether  the  line  joining  the  inner  cusps  of  the  slow 
locus  is  necessaryo   The  complete  answer  would  be  given 
by  evaluation  of  the  Green's  function  for  the  hyperbolic 
system  (2)o   The  relevant  parts  of  the  characteristic 
locus  would  either  carry  singularities  of  the  Green's 
function  or  bound  regions  of  different  analytical  formo 


20  - 


i|..   One -Dimensional  Problems 

The  characteristic  form  of  the  one-dimensional  case  of 
the  nonlinear  system  (1)  or  of  the  linear  system  (2)  is  ob- 
tained by  taking  an  appropriate  linear  combination  of  the 
equations  to  obtain  a  single  equation  in  which  ail  quantities, 
p,  u,   and   B,  are  differentiated  along  the  same  direction  in 
the   (x,  t)  plane.   This  direction  must  be  one  of  the  charac- 
teristic directions.   Six  such  characteristic  equations  can 
be  obtained,  one  for  each  characteristic  direction,  to 
replace  the  six  original  equations.   For  the  linear  system 
with   B   constant,  these  characteristic  differential  equations 
can  be  Integrated  to  yield  six  Rlemann  invariants,  thus  obtain- 
ing an  explicit  solution  for  the  general  initial  value  problem. 
For  the  nonlinear  system,  integration  of  the  characteristic 
equations  is  possible  only  in  special  cases.   This  can  be  done 
if   B   is  either  parallel  or  perpendicular  to  the  x-axls   (  x 
is  the  significant  space  coordinate  in  this  one-dimensional 
problem).   For   B  arbitrarily  oriented,  the  system  can  be 
reduced  to  the  integration  of  one  ordinary  differential 
equation  in  the  special  case  of  a  simple  wave;  i.e.  with  Initial 
data  selected  such  that  there  is  propagation  along  only  one  of 
the  six  characteristics. 

We  denote  by  u   and   B   the  x-components  of  u  and   B. 
■^    n        n 
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Prom   dlv  B=  0,  we  have   B  =  constanto   We  write   u*   and   B" 
for  the  two-dimensional  vector  components  perpendicular  to 
the  X-axis  0   We  write  the  characteristic  equations  of  (1)  in 
terms  of  the  differentiation  operator, 

(8)  d  =  ^  -  (c  -  u^) 


^t         n'  ^x 

If  neither   B   nor   B'   vanishes,  an  elementary  computation 

yields 

/    B'  X  (B  du   +   cdB)    =0      on     c   =   +  A 
/  n  —     n 

(9)  /  2 

A 
pc^du  +  a^c(-~  =  Ddp  -t^B'CB  du+cdB)  =  0    on  c=  +c^,    +  c^, 
\.  c 

The   transverse   speed    is   given  by     A     =  B   Aip,      and   the   slow 
and   fast   speeds    are   denoted   by      c        and      c^    « 

If     B'    =   0     on  some  x-interval,   we  have     u*    =  0     and 

(10)  pdUj^-  cdp  =   0      ,      c   =  +  a        , 

This  is  the  ordinary  gas  dynamics  result  which  can  be  integrated 
to  yield  the  Riemann  invariant 

(11)  u  -p  /  ■— ^  =  constant  on  c  =  +  a  , 


If   B  =0   on  an  interval,  we  have 
n  ^ 


=  22 


u'   =  constant  on  c  =  0 


(12) 


where 
(13) 


B'/p  =  constant  on  c  =  0 


\   ^n  "^  J    "  =  constant  on  c  =  i  a* 


.2   _   2 


(a*)^   =   a^  +  B^Ap  j 


a*  can  be  expressed  as  a  function  of   p  alone  by  using  the 

Integrated  relation  B'/p  =  constant  on  c  =  0   once  the 

initial  values  of   B  and   p   are  given,—' 

For  the  linear  system,  when  B„   and   B'   are  not  zero, 

*'     '        n        o 


we  obtain 


(14) 


B'  X  (B  u  +  cB)   =  constant  on  c  =  +  A 
o   ^  n  —  n 

PqC^u  +  a^c  (-|-Dp  -  ^B^«  (\^"*'  °B)  =  constant 


on  c  =  +  c  o  ■+  c^ 


instead  of  (9).   Also,  In  the  case   B  =  0,  we  obtain  the 
evident  linearization  of  (12),   However,  in  the  case   B^ =  0, 
the  linear  result  is  entirely  different,  viz. 


This  analogy  with  ordinary  one-dimensional  gas  dynamics 
can  be  greatly  extended.   See  Ao  A.  Blank  and  Harold  Grad, 
Fluid  Dynamical  Analogies.   NYO  -  61+86,  Notes  on  Magneto- 
Hydrodynamics   VII,  Institute  of  Mathematical  Sciences, 
New  York  University  (July  15,  1958). 
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p  U    "CD 

^o  n    ^ 


constant  on  c  =  +  a 


o 


(15) 


B  u«  +  cB'-  constant  on  c  =  +  A 
n  -=  n 


The  second  relation„  on  the  characteristic   c  =  +  A   replaces 

the  nonlinear  result   u' =  B'  =  0„  We  see  that  the  restriction 

B' =  0   In  the  nonlinear  case  is  much  more  stringent  than  the 

linearization  about  B' =  0  which  allows  a  non-zero  perturbation 

o  -^ 

B'o   The  linear  result  (15)  can  be  verified  more  rigorously  by 
taking  an  appropriate  limit j,  B'-=— j»0  ^    in  the  general  nonlinear 
case  (9)o   It  is  instructive  to  examine  this  limit  for  the 
linearized  equations ^ 


^t-    -^  Po  Jx 


n 


=   0 


(16) 


< 


Poat    ^o  ^x 


B 


<)B' 


B 


n  y)X 


au« 


n  /^x 


=   0 


>B' 


^  X 


n 


By  ignoring  the  right  sides  of  these  equations,  we  obtain  two 
decoupled  wave  equations  ^  one  for  p  and  ii   with  a^  as  the 
propagation  speed j,  and  the  other  for  u'  and  B»  with  A^  as 
the  propagation  speedo   If  B'   is  considered  to  be  small,  the 
right  hand  terras  represent  a  small  amount  of  coupling o 
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Transverse  Wa 


1/ 

ves— ' 


We  take  the  linear  system  with  B    constant,  (2\   and 
look  for  special  solutions  satisfying  the  following  restric- 
tions 


< 


P   =   0 
div  u  =   0 


(17) 

u-'B  =   0 
o 

^B»B^  =   0 

Under  these  conditions,  curl  B  x  B^  =  (B  «7)B  and 

curl  (u  X  B  )  =  (B  •V)uo   Taking  the  direction  of   B^  as 


the  X-axis,  we  obtain 


()B       _       It-,  1  du 


This  is  a  one -dimensional  wave  equation  system  representing 
the  undistorted  propagation  of  waves  at  a  speed  A^=  B^/s/iip~ 
in  the  direction  of  B^  |   y   and   z   enter  merely  as 
parameters o   The  problem  can  be  solved  separately  for  each 
Individual  magnetic  line  in  terms  of  the  initial  values  on 
that  line.   In  particular^  any  distrubance  which  is  initially 


The  special  solutions  given  here  were  first  presented  in  a 
series  of  lectures  at  General  Electric  in  Philadelphia^ 
July  1957 o 
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found  within  a  certain  flux  tube  (of   B  )  will  forever  remain 
within  the  same  tubeo   For  this  special  class  of  solutions  j, 
Figo  i;  (without  the  dotted  lines)  represents  the  propagation 
of  the  entire  disturbance j,  not  only  the  wave  front «   It 
follows  by  inspection  that  the  solution  of  (18)  satisfies 
(17)  for  all  time  provided  that  it  does  so  initially o   The 
term  transverse  is  clearly  appropriate  for  these  waves  since 
both  u   and   B  are  perpendicular  to  the  direction  of 
propagationo 

We  have  exhibited  a  class  of  special  solutions  of  the 
system  (2)  which  satisfy  the  special  one -dimensional  wave 
equation  (l8)<,   Next  we  show  that  the  two  variables 


(19) 


1   =  B^«curl  B 
"       o 


Co  =  B  "curl  u 
o 


satisfy  the  same  system  (I8)  if   B  and   u   are  taken  from 
an  arbitrary  solution  of  the  full  system  (2)o   Performing 
the  indicated  operations  on  the  last  two  equations  of  (2), 
we  easily  obtain 


(20) 
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For  any  solution  of  the  full  system,   j   and   (.<■'   propagate 
one-dlmenslonally  at  the  apeed  A   , 

Strangely  enough,  the  special  solutions  which  satisfy 
(18)  are  J  in  a  certain  sense,  more  than  could  have  been 
expected  from  casual  inspection  of  the  original  system  (2)o 
Since  u  and  B   in  (I8)  are  each  two-component  vectors, 
the  characteristic  "cone"  (line,  in  this  case)  is  counted 
twicej   (18)  is  a  fourth  order  system.   The  reason  why  this 
can  occur  will  become  clear  later. 
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60   Compressive  Waves 

We  again  take  the  linear  system  (2)  and  this  time  look 
for  equations  which  involve  only  the  four  quantities  (1?) 
which  were  taken  to  be  zero  in  the  last  sections   We  write 


(21) 


a 

13 

div   u 

<T 

= 

P/Po 

P 

= 

Vl^o 

a2 
0 

= 

B?/^Po 

and  obtain,  directly  from  the  full  set  (2), 


(22) 


-|S  +  a2^f  +  A2ip   = 

Jt       O   ^  X 


0  ,  A 


dji         ^j         hz 


If-- 


^u 

X 


=      0 


It  is  easily  verified  that  this  is  a  fourth  order  system  and 
has  as  its  two  characteristic  cones  the  slow  and  fast  cones 
of  the  complete  system  (2),   A  particularly  symmetric  form  is 
obtained  if  a   and  u   are  eliminated? 
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H^        °  ^^2    ay^     ^z2  °  ax2    ^y2  az2 


Jt2  "  ^°^^y2    ^z2'      °   ax2    ^y2   ^^2 


(23)    < 

,2«      .   ,2^   ,2. 

Note  that  the  second  equation  has  a  two-dimensional  Laplacian 
operating  on  0~ ^ 

The  adjective  "compressive"  will  be  used  to  refer  to  a 
combination  of  slow  and  fast  waves  even  though  we  shall  see 
later  that  there  are  special  cases  (indistinguishable  from 
transverse  waves)  with  div  u  =  Oo 

The  system  (21)  which  is  satisfied  by  the  four  variables 

(p,  div  u,  u  ,  B  )   taken  from  an  arbitrary  solution  of  the 

full  system  (2)  is  analogous  to  (20)  for  the  variables 

(B  ' curl  B.  B  'curl  u)o   Together,  one  can  expect  the  combined 
o       '   o  ^_.     » 

system  (19)  and  (21)  to  be  equivalent  to  the  original  sixth 
order  system;  this  will  be  verified  in  the  next  section. 
Special  solutions  of  the  full  system  which  can  be  interpreted 
as  compressive  type  (analogous  to  transverse  type,  (l8))  are 
obtained  as  solutions  of  the  fourth  order  system  (21)  by 
adjoining  the  conditions 

=   B  »  curl  B  =   0 
o 

(21+) 

U)  =  B  •  curl  u  =   0 
o 

Since   j   and  u)  satisfy  the  system  (20),  they  are  identi- 
cally zero  if  they  are  taken  zero  initially. 
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7 o   Resolution  of  an  Arbitrary  Wave  Into  Comp onen t s 

The  results  of  the  two  previous  sections  can  be  summar' 
Ized  In  the  following  tables 


Variables 

Transverse 

Compressive 

p,  div  u,  u^,  B^ 
(curl  B)^,  (curl  u)^ 

zero 
satisfy  (20) 

satisfy  (22) 
zero 

Table  I 

The  separation  of  an  arbitrary  disturbance  into  the  two  com- 
ponents transverse  and  compressive  has  been  accomplished,  but 
only  by  Introducing  nevj  dependent  variables.   Given  an  arbi- 
trary initial  state   (p,  u,  B),  one  can  compute  the  values  of 
the  auxiliary  variables   (PjOt^u  ,  (3^  j,0J)   and  then  find  the 
time  variation  of  the  two  groups   (p,S^y»  P^  ^^^      (jj^) 
independently 8  We  shall  now  attempt  to  describe  this  sepa- 
ration in  terms  of  the  original  variables.   This  amounts  to 
an  inversion  of  the  relations  (19)  and  (21)  to  find   (p,  u,  B) 
when  given   (p,  a  u„,  B  ,  jyW)  o   This  is  nontrivial  since  it 
involves  inversion  of  differential  operators. 

In  terms  of  the  auxiliary  variables,  the  transverse 
component  is  defined  by   p  =  a=u=(3  =  0  with   j   and 
l5  given.   The  transverse  component  in  the  original  variables 
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is  given  by   p=u   =B=0  with  u'   and   B'   (the  compon- 
ents  perpendicular  to  B  )   to  be  found  from 


div'  u»   =   0  fdiv'  B«   =   0 


(25)   -{ 


curl'u'   =  t'^/lB^I  =given  curl«B»   =  j"/|B^j  =  given 


The  accents  (')  denote  two-dimensional  operators  in  a  plane 
X  =  constant.   These  differential  equations  determine   u' 
and   B'   uniquely  in  each  plane   x  =  constant  if  a  regularity 
condition  (e.g.  boundedness)  is  imposed  at  infinity.   Exactly 
the  same  situation  occurs  for  the  compressive  component  except 
that  the  inhomogeneous  terms  occur  in  reverses 

I  div'  u'  =  a-.^u^/^x  =  given  f  div'  B»  =-i^^/^x  =   given 
(26)  \ 

I  curl'u'  =   0  I  curl'B'  =   0 

Both  the  transverse  and  compressive  components  have  now 
been  uniquely  determined,   P,  u^,  and   B   directly,  and  u' 
and   B'   by  (2^)  and  (26).   Furthermore,  the  sums  of  the  two 
components  of  u'   and  of   B'   are  equal  to  the  originally 
given  values.   For,  the  difference  between  u'   and  the  sum 
of  its  two  components  is  harmonic  in  each  plane   x  =  constant 
(i.e.  both  the  operators   div'   and   curl'   annihilate  It), 
and  a  vector  which  is  harmonic  in  the  whole  plane  and  bounded 
at  infinity  vanishes  identically. 
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We  have  obtained  the  resolution  Into  components  of  the 
original  variables.   Unfortunately,  this  resolution  seems  to 
violate  the  domain  of  dependence  properties  which  must  hold 
for  solutions  of  a  hyperbolic  system.   For  example,  consider 
initial  values   (Ps  u,  B)  which  are  identically  zero  outside 
some  bounded  domain.   The  solution  of  the  original  system  (2) 
will  also  be  identically  zero  outside  a  certain  bounded 
domain  which  grows  with  timeo   The  two  sets  of  variables 
(p,  ttj,  u  J,  B  )  and   (j,,^')  behave  similarly  (but  over  differ- 
ent domains)  since  they  also  satisfy  hyperbolic  systems.   How- 
ever, the  transverse  or  compressive  components  of  u'   and   B' 
are  obtained  by  solving  (25)  and  (26)  (essentially  equivalent 
to  Poisson's  equation)  J,  and  they  will  in  general  be  different 
from  zero  over  the  entire  plane   x  =  constant,  even  though 
this  is  not  true  of  the  auxiliary  variables.   Of  course,  the 
sum  of  the  transverse  and  compressive  components  of  u   and 
B  must  still  be  zero  at  large  enough  distances.   In  particu- 
lars, the  resolution  of  the  initial  values   (p,  u^  B)  has  the 
property  that  the  transverse  and  compressive  components  are 
equal  and  opposite  outside  the  disturbed  domain.   We  have  the 
strange  conclusion  that  the  transverse  and  compressive  compon- 
ents which  Initially  cancel  each  other  must  continue  to  do  so 
for  some  finite  time  until  a  vjave  from  the  Initially  disturbed 
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region  catches  up.   Thus,  these  special  initial  values  (viz, 
two-dimensional  harmonic  vectors  in  each  plane  x  =  constant) 
must  propagate  in  exactly  the  same  way  as  solutions  of  the 
transverse  system  which  has  a  speed  of  propagation  A    in 
the   x-direction  and  as  solutions  of  the  compressive  system, 
with  the  complicated  speeds  of  propagation  described  earlier. 
This  is  easily  verified.   As  a  matter  of  fact,  it  was  noted 
in  (17)  and  (I8)  that  even  more  general  sets  of  data  (  u'  and 
B'   not  necessarily  harmonic)  propagate  one-dimens lonally  at 
the  speed   A  « 

The  fact  that  the  resolution  of  initial  values  of  u 
and   B  must  be  done  in  the  large  is  a  consequence  of  the 
property  that  a  fast  or  slow  wave  front  will,  in  general, 
leave  behind  it  a  transverse  residue  which  then  propagates 
one-dimensionally ,   The  significant  point  of  this  analysis  is 
that  one  can  compute  a  pr 1 or 1  what  this  residue  will  be|   the 
transverse  component  of  the  Jjiltlal  u   and   B  gives  this 
residue  on  each  magnetic  line.   In  general,  this  residue  will 
drop  off  as   l/(y   +  z  )   for  large  distances. 

One  might  be  tempted  to  seek  a  resolution  in  a  bounded 
domain  by  solving  the  equations  {2S)    and  (26)  in  such  a  do- 
main and  applying  boundary  conditions.   However,  in  order  for 
a  set  of  initial  values  to  be  termed  "transverse".  It  is 
necessary  that  the  transverse  jump  conditions  hold  at  the 
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boundary  of  the  disturbed  domaino   For  example,  the  condition 
div  u  =  0  would  have  to  be  supplemented  by   u  -  0^  etc o   One 
can  extend  Table  I  to  obtain  conditions  which  guarantee  that 
a  given  bounded  component  is  either  transverse  or  compresslveo 
However,  It  is  not  possible  to  expand  arbitrary  initial  data 
as  a  sum  of  such  modes.   The  essential  point  is  that   the  dif- 
ference between  the  given  data  and  the  sum  of  its  transverse 
and  compressive  components  is  harmonic,  but  in  a  bounded  domain, 
a  harmonic  vector  need  not  be  zeroo   As  an  example,,  consider 
initial  data  in  a  bounded  domain  with  the  property  that  the 
jumps  at  the  boundary  of  the  domain  are  transverse  but  the 
Internal  initial  data  are  note   The  "dominant"  wave  front  will 
be  a  transverse  one,  but  it  will  be  preceded  by  a  fast  front 
which  carries  a  higher  order  jumpo 

There  is  an  alternative  separation  of  the  general  solu= 
tion  of  the  system  (2)  into  waves  of  transverse  and  compres= 
slve  type  which  is,  in  a  certain  sense,  simpler  than  the 
resolution  Into  modes  just  described,,   We  replace  the  system 
(2)  by  the  compressive  system  (22)  in  the  variables 
(a^  Pj  u  5  B  )  together  with  the  following  exact  equations 
for  <)uV<)t   and  ()BV<^t3   taken  directly  from  (2)g 

Po  If-  -  ^^VV)  B'  =  ^  curl  B^xB^  ^  a^V'p 


(2?) 


^-  -  (B  °V)u'   -   0 


o 


3k 


The  solution  of  a  given  initial  value  problem  is  obtained  in 
two  steps.   First   (a,  p,  u  ,  B  )  are  obtained  from  (22)  using 
the  given  initial  values «   Then  these  variables  are  inserted 
into  (27)  which  is  solved  for  u'   and  B'  using  the  given  in- 
itial values  of  these  variables.   This  is  a  one-dimensional 
wave  equation  for  u'   and  B' ,  but  with  a  forcing  term,.   The 
solution  to  (27)  is  the  sura  of  a  special  solution  which  takes 
homogeneous  initial  values,  u'  =  B' =  0,  and  a  solution  of  the 
homogeneous  equation  with  the  given  initial  values.   The  homo- 
geneous solution  will  have  the  one -dimensional  domain  of 
dependence  of  a  transverse  wave.   The  special  solution  will 
have  the  domain  of  dependence  of  a  compressive  wave  since  this 
is  the  region  in  which  the  forcing  terms  differ  from  zero. 

The  domain  of  dependence  complications  associated  with 
the  previous  resolution  into  modes  have  disappeared.   On  the 
other  hand,  although  the  procedure  just  described  gives  a 
succinct  description  of  the  wave  propagation  features  of  the 
system  under  study,  it  does  not  offer  a  separation  into  modes o 
Specifically,  the  separation  at  t  =  0  does  not  agree  with  the 
separation  that  would  be  obtained  using  the  same  procedure  at 
a  later  time.   To  obtain  the  actual  transverse  and  compressive 
components,  one  must  split  the  initial  values  of  u'  and  B'   in 
the  complex  way  described  above  rather  than  take  them  zero  for 
one  of  the  components.   It  is  clear,  from  this  analysis,  that 
any  splitting  of  initial  values  will  leave  unaltered  the  pro- 
perty of  the  sum  that  it  is  zero  outside  the  compressive  domain 

of  dependence. 
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8o   The  Limiting  Cases   A  »  a   and   A  «  a. 

If  the  Alfven  speed  Is  either  very  large  or  very  small 
compared  to  the  gas  sound  speed,  there  Is  a  fiirther  separation 
of  solutions  possible,  vlZc  between  the  slow  and  the  fast 
modes. 

First  we  examine  the  normal  speed  and  characteristic 
loci  to  obtain  the  limiting  behavior  of  the  wave  fronts 
(which  always  separate,  for  any  relation  between  A   and   a). 
Prom  the  characteristic  determinant,  (I4.),  we  find  the  slow 
and  fast  roots 

2     2    2    2 
A  »  a  :   c   =  A  ,   a  cos  0 

A  «  a  s   c   =  a  ,   A  cos  0 

where   0   is  the  angle  between  B  and  the  normal,   n.   In 
each  of  the  two  limiting  cases,  the  fast  normal  speed  locus 
is  a  sphere  centered  at  the  origin,  and  the  slow  locus 
consists  of  two  spheres  through  the  origin.   Prom  this<,  we 
conclude  that  the  fast  characteristic  locus  is  spherical, 
representing  isotropic  three-dimensional  propagation,  whereas 
the  slow  locus  consists  of  two  points  and  the  segment  joining 
them,  representing  one -dimensional  propagation  in  the 
direction  Bo 

The  complete  description,  including  the  transverse  wave, 
is  as  follows,   Por  A  »  a,   there  is  an  isotropic,  three- 
dimensional  spherical  wave  (fast)  with  A  as  the  propagation 
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speed J    and   there    are   two    one -dimensional  waves,    one    (trans- 
verse)   at   the    speed      A      and    the    other    (slow)    at    the    speed      a. 
For     A  «   a,      there    is    a  three-dimensional   spherical  wave 
(fast)    at    the    speed      a,       and    two    superposed   one -dimensional 
waves    (slow   and   transverse)    at   the   speed     A. 

A  slightly  finer  examination  shows  that  the  cusped  part 
of  the  slow  characteristic  locus  shrinks  to  a  point  as  A 
becomes  either  large  or  small  compared  to  a.  This  means 
that  it  takes  a  relatively  long  time  for  the  slow  mode  to 
deviate  from  one-dimensional  propagation.  Of  course,  the 
approximations  made  here  cannot  be  expected  to  be  valid 
uniformly   for    long   times. 

In  order   to    investigate    a  possible    separation  between 
slow  and  fast  wave  motions    (as   distinguished  from  wave 
fronts),    it    is  necessary  to   turn  to   the   differential  equa- 
tions   {22)    for   the    compressive    or   slow  plus  fast    component. 
Note   that      p     and     B        have   been  made   dimenslonless     using 

p        and     B        as    scale   factors,    (21),    but      u^     has  been   left 
'^o  o  '  '  x 

dimensional   since    to    choose    either      a        or     A        as   the    scale 

o  o 

factor  might  prejudge  the  result. 

A  complete  separation  of  initial  conditions  can  be 

obtained  by  Fourier  transform,  expanding  in  powers  of   a  /A 

or   A  /a  .   In  the  case   A  >>  a  ,  the  relative  magnitudes 
0^0  00*  ^ 

of   o~,  p,  and  u   are  independent  of  the  wave  number,   k, 

to  the  lowest  order  in   a  /A  j   the  separation  of  modes  can 

o'  o'        ^ 
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therefore  be  performed  locally  at  each  point  In  space.   For 

A  «  a  ,   the  separation  into  modes  involves  k  even  to 
o     o 

the  lowest  order  in  A  /as   consequently  the  separation 
can  be  made  only  in  the  large o   One  must  solve  differential 
equations  in  order  to  separate  the  initial  values  into  slow 
and  fast  components.  Just  as  was  the  case  for  some  of  the 
variables  in  the  separation  between  transverse  and  compres- 
sive modes o 

These  results  can  be  obtained  more  quickly  by  heuristic 
arguments  stemming  directly  from  inspection  of  the  differ- 
ential equations.   Specifically,  the  fast  mode  is  obtained 
by  formally  setting  the  slow  sound  speed  equal  to  zero  and 
the  slow  mode  by  letting  the  fast  speed  approach  infinity. 
In  essence,  what  we  do  is  present  the  limiting  differential 
equations  whose  solution  amounts  to  inversion  of  the  approxi- 
mate Fourier  transform  obtained  by  taking  a  /A   or  A  /a 

•^      °   o   o       o   o 

small o 

Consider  the  case   A  »  a  .   Setting   a  =  0   in  (22) 

o  o  ^        o 

we    obtain,    for   the    fast   mode 


3^*°=° 


<  7t       ~o 

(29) 


4#  +  A^AP   =   0 


t~=   ° 
Jt   ""   ^  ■  ^  =   ° 
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We  interpret   <)u  / ^\.   =  0   as   u  =  0   since  this  implies 


that  the  value  of   u   propagates  slowly  (with  zero  speed, 
to  this  approximation).   By  the  same  argijment,  we  conclude 
that   p  =  C  from   -a  =  d^/h^   =  ^cr/^t.      The  system  (29) 
can  now  be  put  in  the  form 


(30) 


A  »  a 
o     o 


fast  mode 


For   the    slow  mode  we   let     A 


oo      in    (22)    and   obtain 


(31) 


if--=° 


P   =   0 
TT—  +  a     V"  =  0 


c)u 


V    ^--^  =  0 


¥e  Interpret   Ap  =  0   as   p  =  0   in  the  infinite  domain 
and  obtain  the  simplified  system 


is:  + 


<^u. 


(32) 


P  =  0 

()u 

a  =  -f — 
^x 


X 


=  0 
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A  »  a 
o     o 


sloxif  mode 


We    see   that      u        is    carried  by    the    slow  wave    alone,    and      p 
is   carried  by   the    fast  wave;      cT    and     a      are    carried  by  both. 
By    inspection  of    (30)    and    (32)   we    find   the    following  resolu- 
tion of   initial  values: 

s   _ 
u     =  u 

X  X 


(33) 


P'  =  0  [A^  »   a^] 


C®   =   cr  -    p 


j£_^^x 

It       7^ 


-   a 


The  initial  value  of   <)p  /)t   is  required  in  addition  to  the 

f 
value  of   p   in  order  to  solve  the  second  order  wave  equa- 
tion for   p-^. 

Summarizing,  we  see  that  the  modes  separate  locallyo 
The  longitudinal  component  of  the  magnetic  field  pertxxrba- 
tion  is  carried  Isotropically  at  the  Alfven  speed  and  the 
longitudinal  component  of  velocity  is  carried  one-dimension- 
ally  at  the  gas  speed;  the  condensation  <r,      is  carried  by 
both  waves.   The  slow  wave  is,  in  some  sense,  gas  dynamical, 
but  it  is  enormously  modified  by  the  presence  of  the  magnetic 
field. 

It  is  easily  seen  (cf.  equation  (26))  that  the  remaining 
components  of  vleocity  and  magnetic  field,   u'   and  B',   are 
carried  by  the  fast  wave.   In  other  words,  the  transverse 
residue  is  left  by  the  fast  wave  alone. 
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Now  we  turn  to  the  more  difficult  case,   A  «  a  . 

'^         o  o 

It    is    convenient   to   rewrite    the    system    (22)    in  the 
following  form ; 


+  a  =  0 


(3i^) 


i-f  -   a^  A  cr    =  A^  A  p 
Jt^  ° 

4i   „  A^  tf  =  A^  -^    (^+  ^) 
dt  o   dx  o   ,^x      ^    2       ^^2 


where 


(35) 


[0  = 


*  = 


^P 


AU       -    V- 

X        dx 


It  is  convenient  to  consider  the  terms  on  the  right  in  (3^) 

as   small   coupling   terms |    these    terms   are   removed  when  we 

examine    the    fast  wave      (A     =  O)      or   the    slow  wave      (a      -— >   oo 

o  o 

which  implies   ^^^=0   or  cr    =  O)  o      To  this  approximation 

(and  in  these  variables)  the  system  splits  into  fast  and  slow 
modes; 


(36) 


)t^ 


a  ACT  =  0 
o 


0  =  \Jr  =  0 


A  «  a 
o     o 


fast  mode 
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(37) 


^  =  0 


-  a2  ^  =  0 
o  dn 


A     «  a 
o  o 


slow  mode 


The    complication  occiirs    (as   in  the    separation  of   transverse 
from  compressive  waves    -   of.    section  7)   when  we   return  to    the 
original   variables   via    (35) o      For  the    fast  mode,      p      and     u 

■A. 

are  obtained  by  solving  the  Poisson  equations 


-  d_o'  .  d   o" 

2      — 2" 

<)y    dz 


(38) 


_  da 
^u  -  ■^— 


For   the    slow  mode,      p      and     u        are   obtained  from 


Ap    =    ^ 


(39) 


^u^  =  ^ 


The  given  values  of   Ap   and   au   are  split  into  a  fast 
component  as  in  (38)  with  the  remainder  as  the  slow  component, 
The  fast  components  of   Ap   and   Au^  propagate  three- 
dimenslonally  and  the  slow  component,  one - dime nslonally. 

It  is  not  obvious  (and  is,  as  a  matter  of  fact,  true 
only  to  a  certain  approximation)  that  the  sum  of  the  sepa- 
rated modes  cancels  out  sufficiently  far  from  the  disturbed 
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region.      In  order   to    verify  this,    it    is   convenient   to 
consider   an  alternative   separation  into   waves   of   approxi- 
mately  slow   and  fast    type   similar   to   the   one   described  at 
the   end   of   the   last    section.      We    solve    the    system    (36)    and 
(37)    in  two   steps.      First      a     and      cr     are   found  from   (36). 
Then    (3?)    is    solved   for    the   variables       AP      and      Au 
(rather  than  for     0     and     »|/ )      after   inserting   the    known 
values    of      a      and      Ol      The    solution  of   this   inhoraogeneous 
system   can  be   described  as  the    sum   of    a  homogeneous    solution 
which   takes   the   Initial   values    (propagating   one-dimensionally ) 
and   a   special    solution  with  zero  initial   values    (propagating 
three-dimensionally   as   does    the    forcing   term)c      This    is   not 
the    separation  into   fast   and   slow  modes.      However,    with  this 
splitting,    it    is   easy  to    see   that       Au        and      ^p      are,  for 
all    time,    the    Laplacians    of   fiinctions   which   vanish   outside 
the    domain   of   dependence.      For    the   homogeneous    solution, 
this  property    obtains   initially   and  will    continue    to   do  so 
as   the    solution  of   a   one-dimensional  wave  equation    (cf.    the 
argument   in  section  7)0      For   the   special    solution,    the 
conclusion  follows  provided  that    the    inhomogeneous   term  is 
itself  the  Laplacian  of  a   function  which  vanishes   outside 

the   domain  of   dependence.      These    inhomogeneous   terras   are 
2  2 

(^  l()      C  ^    r) Cfx  ()  /c)ai     _ 

Jtll       o3^^y2        ^^2  H^o  o  ^y2       j^2j 
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To  the  approximation  that  A  «  a  ,   the  statement  is 
^^  o     o* 

verified.   It  is  possible  to  take  a  slightly  different 
approximation  for  the  slow  mode  than  (37),  viz.  the  (exact) 
last  two  equations  of  (3^)>  in  this  case  the  domain  of 
dependence  cancellation  would  hold  exactly,  but  the  separa- 
tion of  modes  would  only  be  approximate. 
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9o   Applications 

Perhaps  the  two  most  Interesting  feattires  of  the  wave 
solutions  discussed  here  are  the  complex  wave  front  patterns 
which  can  arise,  especially  from  the  slow  wave,  and  the  fact 
that  one-dimensional  propagation  can  occur  in  a  three-dimen- 
sional problem.   In  this  section  we  discuss,  very  briefly, 
some  possibilities  of  observing  these  effects. 

In  a  liquid  metal  one  will  have  A  «  a   with  the  attain- 
able magnetic  fields.   The  magnetic  field  and  velocity  perturba- 
tions carried  by  the  fast  wave  are  quite  small,  on  the  order  of 
the  condensation,  cr  =   p/p  .   For  the  same  reason,  the  slow 
residue  left  by  the  fast  wave  will  be  small.   Consequently, 
one  should  expect  to  observe  coincident  slow  and  transverse 
modes,  propagating  as  one-dimensional  disturbances.   The 
finite  conductivity  will  exert  a  very  large  influence  unless 
the  apparatus  is  large  and  the  magnetic  field  very  strongo 
To  observe  the  transverse  or  slow  waves  at  all,  the  product 
of  the  magnetic  field  in  gauss  and  the  apparatus  size  in 
centimeters  should  exceed  about  10   for  mercury  or  2  x  10^ 
for  liquid  sodixom.   To  detect  the  cusps  developing  from  the 
slow  wave  is  much  more  difficult |  for  this,  the  apparatus  size 
is  not  so  ci^itical,   but  the  magnetic  field  should  be  at 
least  one  million  gauss  in  mercury  or  10"^  gauss  in  sodium. 

The  one-dimensional  propagation  is  of  particular  interest 
because  it  represents  a  transport  of  energy  which  is  not 
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radially  attenuated.   Probably  the  most  interesting  example 
is  the  propagation  of  signals  in  the  geomagnetic  field.   The 
theory  presented  here  offers  only  a  tentative  indication  of 
what  may  be  expected.   In  addition  to  the  poor  conductivity 
and  low  degree  of  ionization,  there  are  important  effects 
associated  with  finite  mean  free  path,  finite  gyro  radius, 
and  possibly  finite  Debye  radius—'  ,   For  most  of  the  range 
of  altitudes,  the  approximation  A  »  a   is  valid.   There 
exist  two  one-dimensional  modes  propagating  along  magnetic 
lines;  one  at  the  speed  a   and  the  other  at  the  speed  A  » 
The  question  of  general  wave  propagation  into  non-uniform 
magnetic  fields   B   has  not  been  treated  here.   However, 
the  propagation  of  wave  fronts  has  been  completely  ana- 
lyzed—^ and  is  in  this  case  strictly  one-dimensional. 
There  is,  of  course,  a  distortion  of  the  wave  front  caused 
by  the  varying  speed  of  propagation  (see  Figure  6).   It  is 
easy  to  compute  that  this  distortion  will  not  cause  the 
amplitude  of  the  signal  to  vary  appreciably  (i.e.  by  more 
than  one  order  of  magnitude)  in  propagating  over  a  distance 
comparable  to  the  earth's  size. 

A  further  question  which  must  be  answered  is  the  transfer 

Some  of  these  effects  have  been  taken  into  account  to  a 
rough  approximation  in  an  analysis  of  "whistlers"  by 
L.  R.  0."  Storey,  Trans.  Roy.  Soc,  2jj^,  113  (1953). 

2 

See  footnote  on  page  6. 
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of  these  magnetohydrodynamlc  signals  in  the  ionosphere 
through  the  atmosphere  to  the  earth.   The  three  modes 
existing  in  the  conducting  medium  are  coupled  to  two 
which  exist  in  air;  viz.,  ordinary  sound  waves  and  vacuum 
electromagnetic  waves.   The  reflection  and  transmission 
coefficients  across  a  discontinuity  siorface  can  be  easily 
obtained  (this  has  not  yet  been  done),  but  the  problem  of 
a  layer  of  variable  properties  will  present  more  difficulty. 
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It  is  worth  noting  that  the  electromagnetic  wave  will  be 
observed  as  a  quasistatlc  magnetic  field  rather  than  as  a 
conventional  propagating  wave.   This  is  a  consequence  of  the 
fact  that  the  speed  of  light  is  much  greater  than  the  propaga- 
tion speeds  in  the  conducting  medium.   The  electromagnetic 
disturbance  in  the  atmosphere  can  be  computed  at  each  instant 
as  a  vacuiom  magnetic  field  which  attains  the  instantaneous, 
value  of  the  normal  component  of  B   at  the  interface.   The 
perturbation  of  B   that  is  observed  will  be  transverse  and 
will  be  associated  with  the  speed  A. 

In  a  laboratory  plasma,  the  limitations  of  this  theory- 
are  equally  strong.   However,  it  is  at  least  worth  remarking 
that  the  usual  investigations  which  are  limited  to  the  propa- 
gation of  plane  waves  only  may  very  well  mask  unusual  and 
important  effects. 
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